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The Bogoliubov approximation and the Gross-Pitaevskii equation characterize the effect 
of repulsive interactions on a dilute ideal Bose-Einstein gas in a spherical harmonic trap. 
For large N, the interactions expand the condensate relative to an ideal Bose gas; both 
the speed of sound and critical angular velocity Qci for creation of a quantized vortex 
depend crucially on the interparticle repulsion through the coherence length ^. 



1 Ideal Bose gas 

An ideal Bose gas provides a valuable introduction to the physics of a real dilute 
Bose gas, for it emphasizes the fundamental role of the coherent condensate that 
contains a macroscopic number of particles in a single quantum state. The standard 
example treats noninteracting bosons in a volume V with periodic boundary con- 
ditions and uniform density n = N/V. In the classical limit, the thermal de Broglie 
wavelength Xt = (2TTh'^ /mkBT)^^^ is much shorter than the interparticle spacing 

1 ~ n~^^^. As the temperature falls, however, the thermal de Broglie wavelength 
grows, and the system eventually becomes degenerate when \t ^ I. Equivalently, 
Bose-Einstein condensation in a uniform ideal gas occurs at Tc ~ h^n?/^ /kBm. 

The situation is somewhat different for an ideal Bose gas in a harmonic trap 
(taken as isotropic for simplicity) with Vtrap('") = ^mujQr'^; the corresponding oscil- 
lator length do = {h/muJoy^^ characterizes the size of the (Gaussian) ground state. 
In the classical limit {ksT ^ ^11^0)7 the density follows the Boltzmann distribution 
riciir) oc exp[—Vtrapir)/kBT], which can be rewritten as exp(— r^/2i?|,), with 
Rt = dQ{kBT /IiujqY/'^ the classical thermal radius of the trapped gas (note that 
Rt > do). 

To estimate the temperature Tc for the onset of Bose-Einstein condensation in 
a trap, use the same expression with n ^ N/Rj., so that fcsTc ~ 'h? N'^^^ /mR^ ^ 
U^ojIN'^/'^ /kBTc] equivalently, fc^Tc ^ N^/^htoo > hujo- Since do < Rt, the con- 
densate forms a narrow spike of width do superimposed on the smooth background 
of width Rt, as is seen clearly in the original experiment with a few thousand *^Rb 
atomsid More recent experiments have condensed millions of atoms, with tynical 
parameters iV « 5 x 10^, tJo/27r « 120 Hz, and do ~ 1.9 /im for sodium atomsH 

2 Effect of Repulsive Interactions 

The interparticle spacing is typically large compared to the range of the atomic 
interactions, so that the interparticle potential can be taken as U{r) « C/o^(r), 
with Uo — ^T^afi' jm and a the s-wave scattering length (a k, 2.75 nm for sodium 
atoms). The Bogoliubov approximationB assumes that nearly all particles are in the 
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condensate; at zero temperature, the system forms a dilute Bose gas with na^ <C 1. 
2.1 Quasiparticles and the Coherence Length 

A uniform Bose gas at T = K has only two characteristic single-particle energies: 
the kinetic energy = h^k'^/2m and the "Hartree" energy 

VH{r) = / U{r - r')n(r') = Uqu = inah^n/m. (1) 

Bogoliubov showed that a quasiparticle with wave vector k has an energyil 

p - /oT/ ^ I ^2 _ / ^sk, for fc ^ (phonons), , . 

^ - y^vnek + Efe ~ I ^ n2fc2/2m, for A: ^ oo (free particles), ^ > 

where = Uon/m = ATrah^n/rn^ is the squared speed of sound. The quasiparticle 
is a linear combination of a particle and a hole. For fc — > 0, it is an equal particle-hole 
admixture, but it becomes a pure particle as fc ^ oo, with the cross-over occurring 
at k ~ where ^ = (STran)"^/^ defines the "coherence" length. Note that the 
speed of sound can be rewritten as s = K/27r-\/2f, where n — h/m is the quantum 
of circulation. When na^ ^ 1, it is easy to verify that £, ^ I ^ a. 



2.2 Radius of Interacting Condensate in a Trap 

The presence of the tran^introduces a third characteristic energy, and the Gross- 
Pitaevskii (GP) equatioiJj'B for the nonuniform condensate wave function ^'(r) has 
the form 

(r + ytrap+T///)*-A^«', (3) 

where T — —H^V^ /2m is the kinetic energy operator, Vtrap(r) = ^hLjQ{r / do)"^ is 
the trap potential energy, Vh (r) = Uon{r) = 'iTiah^n{v)/m is the Hartree potential 
energy of one particle with the remaining particles, and is the chemical potential. 

For very low density, the coherence length ^ exceeds the trap size doj a-nd the 
system acts like an ideal Bose gas. As N increases, however, ^ shrinks; when ^ 
becomes comparable with d^, the repulsive interactions begin to predominate, and 
the self-consistent radius i?o of the condensate expands beyond do . 

To estimate the actual radius, note that the kinetic energy per particle is of 
order T ~ h^/mRQ '--^ htoo/R'^, where R = Ro/do is the dimensionless radius of 
the condensate. Similarly, the order of magnitude of the remaining single-particle 
energies are t4rap ^ hujoR^, and 

h^a h?a N hujQ Na , . 

(4) 

where the dimensionless parameter rj = Na/dg characterizes the strength of the 
interactions. For 77 <C 1, the system is effectively ideal, but for 77 ^ 1, the in- 
teractions become crucial (nevertheless, the system remains dilute with a <C I for 
all practical trapped condensates). If (• • •) = J dV'i'* • • • v]/ denotes a condensate 
ground-state expectation value, the total energy E = {T + Vtrap + ^^h) is of order 
~ NhuJo{R^^ + R"^ + rj R^^), and R is determined by minimizing E. 



2 



For small rj, the minimum energy occurs at i? « 1, and the condensate radius 
is just do, with /i « ^fii^o (the ground-state energy of the isotropic oscillator). 
For large 77, in contrast, the kinetic energy is negligible, and the minimum total 
energy occurs for ^ rj. In this limit, a detailed calculationO shows that /i = 
^huiQ and that R = Ro/do — (15r/)^/^, which quantifies the expansion of the 
condensate relative to do. For parameters appropriate to the MIT experimentsp 
the dimensionless condensate radius is i? w 10.2, so that Rq ~ 19.3 fim. 

The resulting "Thomas-Fermi" (TF) approximation neglects the kinetic energy 
entirely, and the GP Eq. (||) then determines the condensate density by the condition 
Vtrap(r)+VH'(r) = fi, showing that the resulting density profile is parabohc. The cor- 
responding central density n(0) = R^/^iradQ = (IS/Stt) N/Rq serves to define both 
the speed of sound = 47ra?i^n(0)/m^ and the coherence length = l/87ran(0) 
for a trapped condensate; as a corollary, the relation ^Rq — dg implies that do is 
the geometric mean of f and Rq- The previous expression s ~ K/2Tr\/2^ then shows 
that the speed of sound in a large trapped condensate increases linearly with Rq, 
reflecting the increased density. For the MIT experimentsp these relations yield 
n(0) « 4 X 10^0 m-3, s « 1.04 cm/s, and ^ ~ 0.187 /zm. 

In the Bogoliubav approximation, most of the particles remain in the conden- 
sate, which requiresB [n{0)a^]-^/^ <C 1. For a trapped condensate at T = K, this 
condition of small total noncondensate number holds for R^a/d^ <C 1. When this 
dimensionless ratio becomes comparable with 1, however, the repulsive interactions 
are so strong that the total noncondensate number becomes of order iV (this situa- 
tion occurs in liquid helium, even at T = K) . Comparison with the TF expression 
rj = Na/do R^ shows that the Bogoliubov approximation fails for a trapped con- 
densate when N ^ (do/a)^, which is far larger than any current experimental value 
(for the MIT trap with sodium atoms, this limit merely requires that N ^ 10^^). 

3 Critical Linear and Angular Velocities 

Landau's original explanation of superfluidity involved the critical velocity Vc] it is 
the speed at which a moving macroscopic impurity can create quasiparticles and 
thus lose energy. A simple analysis shows that Vc is the minimum value of the ratio 
ujk/k, where Wk is the dispersion relation for the quasiparticles. From this perspec- 
tive, a uniform ideal Bose gas has zero critical velocity, for the dispersion relation is 
simply Tik"^ /2m. In contrast, Eq. (||) shows that a uniform dilute interacting Bose 
gas indeed has a nonzero critical velocity, with Vc — s a^/^, arising from the 
presence of the repulsive interactions. 

3.1 Critical Angular Velocities in a Trapped Condensate 

For a dilute trapped condensate, it is natural to take the speed of sound s as the 
critical velocity. It is impractical to shoot an impurity through the condensate, but 
the equatorial speed of a rotating condensate can serve to define a corresponding 
critical angular velocity — s / Rq — h/mdQy/2 — ujQ/y/2. Thus a large vortex- free 
condensate rotating at an angular velocity < Qc should remain superfiuid (this 
expression also indicates that the frequencies of low-lying compressional modes in 
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a trapped condensate are independent of the radius and of order cjq)- 

The low-lying hydrodynamic modes of a large spherical condensate have the 
dispersion relation iVni = uJo[l + n{2n + 21 + 3)]^/^, where n is the^radial quantum 
number and I is the orbital-angular-momentum quantum numberQ For fixed Z, the 
radial wavenumber k is ~ n/Ro, and a corresponding more precise Landau critical 
angular velocity is V^uiq (this value occurs as fc — > oo). 

It is helpful to review the effect of rotation fl — flz on a, sample of liquid 
helium. If the fluid is normal, the microscopically rough walls bring it into solid- 
body rotation Vgb = 17 x r, where r is the distance from the axis of rotation; this 
flow has uniform vorticity, with V x Vsb — 2r2. In the rotating frame, the walls 
are stationary, and the relevant zero-temperature thermodynamic function is the 
"free energy" F = E — ■ h where E is the total ground-state energy and L is the 
total ground-state angular momentum. For a circular cylinder of radius i?o, the free 
energy of a state with one vortex on the symmetry axis becomes lower than that 
with no vortex at a lower critical angular velocity fid = (/t/27ri?§) ln(it!o/C)) where 
^ represents the vortex core radius is a few atomic diameters for superfluid ''He). 

Assuming that a.^imilar expression holds for the creation of a vortex in a large 
trapped condensatePu the lower critical angular velocity is smaller than the trap 
frequency ujo = K/2TrdQ by a factor of order (^/i?o) ln(-Ro/C) ^ 1- addition 
to the compressional modes with frequencies of order ujq, the presence of a vortex 
line introduces new dynamical degrees of freedom associated with "vortex waves." 
At long wavelengths (fc^ <C 1), the classical vortex- wave dispersion relation « 
(Kfc^/47r) ln(l/fc^) immediately suggests low-lying normal modes with k ~ Rq^ and 
frequencies of order Od, which might serve to signal their presence. 

3.2 Analogy with Type-II Superconductors 

These results for flc and Qci are very similar to the thermodynamic critical field 
He = ^o/2'kV2X£_ and lower critical field Hd ~ ($o/27rA^) ln(A/^), where $o = 
h/2e is the flux quantum, A is the penetration length, and ^ is vortex core radius 
(also the superconducting coherence length)u A uniform bulk superconductor be- 
comes unstable with respect to normal metal at the field He, and the formation of 
quantized flux lines (vortices) becomes favorable at the field Hd ■ In a type-II super- 
conductor (one with A > ^/\/2), the thermodynamic instability at He is preempted 
by vortex formation at Hd < He- 

The electromagnetic currents in a charged superfluid cut off the logarithmic 
intervortex interaction potential, screening it exponentially beyond the penetration 
length A. This quantity diverges as the charge on each particle tends to zero, and 
the corresponding "screening" length in a neutral superfluid becomes either the 
radius of the container or the intervortex separation, whichever is smaller. 

In addition, a type-II superconductor ultimately becomes normal at the upper 
critical field Hc2 = <5'o/27r^^, roughly when the vortex cores overlap. Unfortunately, 
the corresponding flc2 = K/27r^^ is unattainably large in superfluid "^He, and me- 
chanical instability for Q > luq may also render it unobservable in a rotating dilute 
trapped condensate (in this trapped condensate in equilibrium could contain 

only relatively few vortices). 
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3.3 Non-Axisymmetric Rotating Traps 

The magnetic field that confines a dilute atomic Bose condensate acts simply as a 
potential Vtrap (r) ; this situation differs greatly from the microscopically rough walls 
of a container for superfluid helium. Thus, a "rotating trap" is meaningful only to 
the extent that it is nonaxisymmetric, for the rotating time-dependent potential 
pushes the condensate, setting it into motion. 

In contrast to the case of an axisymmetric potential such as a circular cylinder, 
the free energy E — D,L of an irrotational vortex-free state in a nonaxisymmetric 
trap decreases with increasing (like — where / is an effective moment of 

inertia) . In the limit of large distortion, / can approach that of solid-body rotation. 
The negative free energy of the vortex-free state typically delays the onset of vortex 
formation, for the vorticity localized in the vortex cores becomes less necessary 
to "mimic" the uniform vorticity of Vgb- This effect is readily verified in simple 
cases.jJ'or a uniform superfluid in a long rotating elliptic cylinder with semiaxes a 
and the preceding expression Od = (K/27r&^) ln(&/^) applies if b — a, but the 
corresponding lower critical angular velocity fici « (K/47r&^) ln(6/f) for 6 ^ a can 
become significantly larger. 

The GP Eq. (^) provides a basis for analyzing a rotating trap, where the con- 
fining potential Vtrap(r) is stationary in the rotating frame (r now denotes the 
coordinate in the rotating frame). The free energy is given by — flLz, where E is 
the ground-state energy (assuming T = K for simplicity). The condensate wave 
function 5* — |4'|e*'^ can be expressed as a magnitude I^Pl = n}^"^ and a phase S 
that determines the superfluid velocity v — {'h/m)WS. The TF limit of a large 
condensate neglects the spatial variation of n, and the free energy becomes 

F sa J(fr {^mnv^ + nVtrap + lUon^ - mnVlz • r x v) . (5) 

In a singly connected container, v can be written as a sum of two contribu- 
tions: vn oc arising from the moving walls and arising from the vortices with 
circulation k. Correspondingly, F separates into a (vortex-free) contribution En 
associated solely with vq and F^ — J (Pr mn-va ■ -I- ^mnv'^ — mnilz • r x 
that also depends on the angular velocity fi, both explicitly and through vo. The 
critical angular velocity fid for vortex creation occurs when first vanishes. 

The integral for E provides a variational expression in the rotating frame; as 
an approximate trial function for the phase, take S = —{h/m)<^, where $ is the 
classical velocity potential for a uniform irrotational fluid (conventionally defined 
so that V — — V$). This velocity field is automatically irrotational apart from the 
vortex cores and satisfies the condition that its normal derivative match the normal 
velocity of the rotating boundary. In addition, V • v = 0, which is the appropriate 
TF limit of the continuity equation V • (fiv) = 0. The actual trapped-condcnsate 
density differs from that for = only because of the Coriolis and centrifugal forces 
of order mJl x v and mn x (fJ x r), respectively. Each of these is small compared 
to the force of the harmonic trap so long as il ^ in which case it is permissible 
to neglect the change in density (apart from the formation of the small vortex core 
of radius « ^). 

As a result, the free energy _Fh for an irrotational condensate in a rotating trap 
has the same form as that for a uniform incompressible irrotational fluid, but with 
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the particle density taken from the solution of the GP equation. Since the actual 
density varies slowly in the TF limit, differs from the classical expression only 
through factors of order unity, representing various moments of the actual density. 
Similarly, the additional free energy arising from the presence of a vortex line 
also differs from the classical expression for a uniform incompressible fluid only 
by factors of order unity, because the TF density profile cuts off the logarithmic 
divergence in the kinetic energy at a distance of order ^, which caitbe identified with 
the classical vortex radius. Consequently, the classical analysisEEl for uniform fluid 
in a rotating elliptical cylinder provides a qualitative guide to the corresponding 
problem of vortex nucleation in a rotating elliptical trap. 

This correspondence remains valid for a multiply connected container, when the 
velocity field v includes an additional contribution vr arising from the quantized 
circulation F = Jk around the various internal boundaries (here, j is an integer). For 
example, consider an incompressible fluid in an asymmetric annular region between 
two nonconcentric cylinders (outer radius R2 and inner radius Ri < R2), with the 
centers displaced a distance d < R2 — Ri- If the system rotates in equilibrium 
at angular velocity around the symmetry axis of the inner cylinder, the critical 
angular velocity for the creation of vortex-free quantized circulation k around 
the inner boundary can be evaluated for arbitrary d. In the symmetrical limit 
{d 0), the result f7„ = (K/27r)(i?| - i??)-^n(i?2/i?i) is similar to that for a 
cylinder; in the small-gap limit, however, the resulting expression fJ^ ~ k/ A-kRiR2 
has no logarithmic factor. These expressions can serve to estimate the corresponding 
quantities for a trapped toroidal condenaa-te (created, for example, by piercing the 
condensate with an off-axis laser beam)0 The circulation-induced deformation of 
the condensate might serve to detect the presence of a persistent current. 
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